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Abstract 

We present a method to compute the genus expansion of the free energy 
of Hermitian matrix models from the large N expansion of the recurrence 
coefficients of the associated family of orthogonal polynomials. The method is 
based on the Bleher-Its deformation of the model, on its associated integral 
representation of the free energy, and on a method for solving the string 
equation which uses the resolvent of the Lax operator of the underlying 
Toda hierarchy. As a byproduct we obtain an efficient algorithm to compute 
generating functions for the enumeration of labeled /c-maps which does not 
require the explicit expressions of the coefficients of the topological expansion. 
Finally we discuss the regularization of singular one-cut models within this 
approach. 
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1. Introduction 

In this paper we consider the ensemble of random Hermitian matrices 



Z N ( S ) = [ exp (-N V V( Xi , g) ) TT 



(xi — Xj) 2 dxi ■ ■ ■ c\xn, (1) 
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for a given polynomial potential 

V(z,g) = Y,9nZ n (2) 

n=l 

of degree 2p with real coefficients g = (gi, . . . , g 2p ) such that g 2p > (we will 
not make explicit the dependence on g of the functions associated with the 
model © unless there is risk of ambiguity). For more than thirty years [l|, 
0, S 0, 0, @] the asymptotic behavior of the free energy 

FN = ~^\nZ N (3) 

as iV — > oo and its relation to the counting of Feynman graphs have been 
subjects of intensive research. However, rigorous proofs of the existence of 
an asymptotic expansion of Fn in powers of N~ 2 were provided only rather 
more recently by Ercolani and McLaughlin [t| and by Bleher and Its j^]. 
These analyses prove the existence of a genus expansion of the form 

F N (g)-F^J2F (k) (s)N- 2k , (4) 
where F^ stands for the Gaussian free energy, 

under the assumption that there is a path g(t) in the space of coupling pa- 
rameters connecting V(z, g) to the Gaussian potential z 2 in such a way that 
V(z,g(t)) is a regular one-cut model for all t. The functions F™(g) are 
important objects because the coefficients of their Taylor expansions at the 
Gaussian point gj? = 5k,2 are generating functions for the enumeration of la- 
beled fc-maps with vertices involving valences 1, . . . , 2p, where 2p is the num- 
ber of nonvanishing coupling parameters g n . The general aim of the present 
work is the characterization of the structure of these functions -F^'^(g). (Inci- 
dentally, for multi-cut models it has been shown jsl, that in general the free 
energy exhibits an oscillatory behavior as a function of N, and consequently 
topological expansions cannot exist.) 
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The large N asymptotics of the matrix model ([T]) is intimately connected 
with the asymptotics of the recurrence coefficients r n A r and s H! n in the three- 
term recursion relation 

xP n , N (x) = P n+liN (x) + S n , N P n , N (x) + r n>N P n -l,N (x) , (6) 

for the orthogonal polynomials P n ,N(x) = x n + a n -ix n ~ l + ■ ■ ■ with respect 
to the exponential weight 

PkAx)Pi,N(x)e- NV ^dx = 5 k>l h k>N . (7) 



In particular [10|, 111] , in the limit 



n -> oo, N oo, — ->■ T, (8) 

the density of zeros of P n> ^(x) reduces to the eigenvalue density of the matrix 
model Z N (g/T), and if Z N (g/T) is a regular one-cut model then 0, [l2[ the 
recurrence coefficients r n A r and s n> ^ can be expanded in powers of N~ 2 

The methods that exploit this relation with orthogonal polynomials to 
calculate the asymptotics of the free energy essentially consist of three steps: 

1. Determine the free energy in terms of the recurrence coefficients. 

2. Obtain the asymptotic expansion of the recurrence coefficients. 

3. Use 1 and 2 to obtain the asymptotic expansion of the free energy. 

There are alternative methods based on solving the Ward identities for the 



partition function (loop identities) [13] (cf. also |14|) or on formulating the 
matrix model as a conformal field theory [lij . However, the structure of the 
expressions for the F^(g) that these methods provide is less suitable than 
ours to compute generating functions for the enumeration of labeled fc-maps. 

In this paper we present an efficient method to calculate the topological 
expansion of the free energy and to characterize its coefficients. For simplicity 
we restrict our analysis to Hermitian models associated to even potentials 
V(A), where A = z 2 . We now preview how our approach performs the three 
steps and point out the differences with respect to other schemes. 



The classical method of Bessis, Itzykson and Zuber [g, [16|, [17|, [18|, [19|, [2 



21| is based on the following identity (step 1) 



v = -^ln(iV!) - - ln^iv + £ (l - lnr ni7V , (9) 



n=l 
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wherein the expansion of r n A r is substituted (step 2), and the asymptotic 
behavior of F^ is obtained by means of the Euler-Maclaurin summation 
formula (step 3). However, some objections to this approach were raised by 
Ercolani and McLaughlin (cf. subsection 1.5 in (7|) due to the use of the 
asymptotic series of T n ]\f clS db uniform expansion valid even for n = 1 as 
iV — > oo. This objection triggered the interest in alternative strategies for 
step 1 (21, 0, [22}] . For example, Ercolani, McLaughlin and Pierce 22] derived a 
hierarchy of second order differential equations to determine the coefficients 
of the expansion (j3J) from those of the the asymptotic expansion of r n ^. In 
our work we use instead the Bleher-Its integral representation ^ 



F N (g) = F% + 



1 -t 



r N)N (g(t)) (r N -i, N {g(t)) + r N+1>N (g{t))) - 



where g(t) denotes the Bleher-Its deformation [2| 

V(X, g(t)) = (1 - l/t)X + V(X/t, g), 1 < t < 00, 



or explicitly in terms of the coupling parameters, 
g 2 (t) = 1 - - + ^, g 2k (t) 



92k 

t k ■ 



k > 2. 



dt. 
(10) 

(11) 
(12) 



For step 2 the standard methods 0, 0] use recursive equations for the 
coefficients (string equations) to determine expansions of the form 



^(g)~$>(T,g)e 2fc , e = l 



(13) 



fc>0 



Bessis, Itzykson and Zuber [17|, |23| formulated the string equation in terms 
of summations of paths over a certain staircase. Some years later Shi- 
rokura [19|, |20J developed a general method to perform these summations 
and characterize the coefficients of the expansion f[T3"j) in terms of the single 
function 

'2n\ 



^.s>=£( 2 r) 

n=l ^ ' 



ng 2n r ■ 



(14) 



To solve the string equation, in this paper we introduce a generating function 
U u ^n{X) associated to the resolvent of the finite-difference Lax operator L of 
the underlying Toda hierarchy \2A 



LP, 



n+l,N 



+ r n . N P, 



(15) 
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and determine the coefficients of (1131) from the function (1141) as rational func- 
tions of the leading coefficient r . This approach is simpler that Shirokura's 
method, is particularly suitable for symbolic computation, can be applied to 
a generic potential, and permits an easy characterization of the asymptotics 
of the Bleher-Its deformation r n jv(g(^)) of the recurrence coefficient. 

Finally, regarding step 3, the Bleher-Its representation f lTUj) allows us to 
express the coefficients of the topological expansion as 

pro 

F (k \g)= / R k (Z,g)d£, k>l, (16) 



where the integrands Rk{£,, g) are rational functions of £ which can be com- 
puted explicitly in terms of the function ffT4j) . 

The layout of this paper is as follows. In section 2 we briefly review the 
basic facts about matrix models which are required to discuss the Bleher-Its 
deformation and the conditions that ensure the existence of the expansion (J3J). 
Then we apply our results to the quartic model and to the sixtic model of 



Brezin, Marinari and Parisi |25[. Section 3 is devoted to the asymptotic 
expansion (fl3|) of the recurrence coefficient for general models V(X, g) and 
their respective Bleher-Its deformations V^(A, g(t)). In particular, we rederive 
in a much shorter way the expressions for the coefficients r^T, g) and r 2 (T, g) 



found by Shirokura [19|, |20| . Section 4 deals with the asymptotics of the free 
energy. From the Bleher-Its representation ( imp and the expansion of the 
recurrence coefficients we obtain the integral expression ffl6|) . We evaluate 
explicitly the integrals for the up to genus 3 in the general case and find, 
except for a coefficient in the expression of F<- 3 \ the same results found by 



Shirokura [19|, [20[ . We also check that the expressions of these coefficients 
for general two-valence models reduce to those obtained using the Ercolani- 
McLaughlin-Pierce method j22|. In the brief section 5 we discuss how to 
apply our method to compute counting maps functions and present explicit 
calculations for two and three valence models. In section 6 we formulate 
a "triple scaling" method to regularize the free energy expansion of a class 
of singular models (the singular one-cut case) and show how the Painleve I 
hierarchy emerges in our approach. The paper ends with a brief summary. 

2. The Bleher-Its deformation 

To calculate the asymptotic behavior of F N (g) as N — > oo using the 
Bleher-Its formula flTTJj) we need the asymptotics of the deformed recurrence 
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coefficients r n ^(g(t)) where g(t) is the Bleher-Its deformation of g. In this 
section we study the action of this deformation on the space of coupling 
parameters. As relevant examples we analyze the Bleher-Its deformation for 
models associated to quartic potentials and to the sixtic potentials of Brezin, 



Marinari and Parisi 25 



The continuum limit of r n ^(g(t)) depends on the number q of cuts of the 
model V(X, g(t))/T. In appendix A we summarize the method to determine 
the number of cuts of a generic hermitian model. The endpoints of a q- 
cut eigenvalue support J = Uj =1 (<x,-, f3j) are the solutions of the system of 
2q equations f l 1 8 1 [) — (II 8 3 j) but, in general, several such systems of equations 
corresponding to different values of q may have admissible solutions for one 
and the same model. Among these candidate solutions, the correct value of 
q is uniquely determined by the additional set of inequalities f ll84p - fll86p on 
the polynomial h(z) defined by 

= h(z) + Oiz- 1 ) as 2^00, (17) 

Wi(z) 

where W\(z) is the branch of the function 



n(*-ai)(z-A) (is) 

\! i=i 

with asymptotic behavior W\(z) ~ z q as z — > 00. In turn, the polynomial 
h(z) is related to the eigenvalue density p(x) by 

h(x) 

p( x ) = -tH- w i,+( x ) for x e J, (19) 

Z7T1 

where wi !+ (x) denotes the boundary value of W\{z) on J from above. 
We restrict our considerations to even potentials of the form 

p 

V{X,g,) = J292^ j , A = ^ 2 , g 2p >0, (20) 
i=i 

where the coupling constants g = (g 2 , (74, ... , g2 P ) run on a certain region G 
of R p . The phase diagram of the corresponding family of matrix models is 
introduced through the decomposition 

G=(jG qi (21) 

9=1 
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where g G G q if and only if g determines a g-cut regular model (cf. ap- 
pendix A). We will refer to G q as the g-cut phase of the family (|2"U|) of 
Hermitian models. For even potentials the eigenvalue support J is sym- 
metric with respect to the origin, and in the one-cut case the endpoints of 
J = (—a, a) are determined by the single equation (11831) . which in terms of 
V reduces to 

5s v Wr^s = 1 - (22) 

z7Ti V A — a z 

Regarding the behavior of a particular model g G G with respect to its 
Bleher-Its deformation, we will consider two cases in our analysis: the regular 
one-cut case, in which g(£) G Gi for all t > 1, and the singular one-cut case 
in which g(t) G G*i for t > 1 but g = g(l) determines a singular model 
(cf. appendix A). 

2.1. The quartic model 
The quartic model 

V(X,g)=g 2 X + gi X 2 (23) 

in the region 

G = {g = (92, 9i) GR 2 : 24 >0}, (24) 

only exhibits g = 1 and q = 2 phases [2a |. The g = 1 phase G\ can be written 
as the union 

G x = Gf } U Gf } , (25) 

where 

= {fe,34)6R 2 :<?2>0,< ? 4>0}, (26) 
G? = {(g2,94)eR 2 :g 2 <0,g A >0,g 2 >-2^g- i }, (27) 

and the q = 2 phase G2 is given by 

G 2 = {(32,34) G R 2 : g 4 > 0,92 < -2^}. (28) 

The phase diagram features the critical curve 

92 = -2y/gl, (29) 

which demarcates the transition line between the two phases. 
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Figure 1: Deformation paths g(t) for the quartic potential. From left to right: a defor- 
mation of a two-cut regular model, a deformation in the singular one-cut case and two 
deformations in the regular one-cut case. 

Consider now the Bleher-Its deformation (cf. figure [T|) 

g(f) = fat)Mt)) = ( t ~\ +92 ^) • ( 3 °) 

If g G Gp then g 2 > and g 2 {t) > for alH > 1. Therefore g(i) G Gi for 
all i > 1. If g G Gf \ then # 2 > -2y/g~4 and 

= > = -2v^)- (31) 

Hence if g G G\ we also have that g(£) G Gi for all £ > 1. On the other 
hand, it is elementary to see that if g G G 2 or is on the critical curve (129]) 
then g(to) is on the critical curve for t = 1 — g 2 — l^fgi- Summing up, 

1. If g G G x then g(i) G G x for all t > 1. 

2. If g G G2 then g(t) crosses the critical curve at to = 1 — 92 — ^■sfgi- 

3. If g is on the critical curve f )29|) then g(t) G Gi for alH > 1. 
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2.2. The Brezin-Marinari-Parisi model 



In [25j Brezin, Marinari and Parisi considered the potentials V(z)/T with 

V(z) = 90z 2 - 15z* + z 6 (32) 

to generate a non-perturbative ambiguity-free solution of a string model. The 
models V(z)/T define a path in the space of coupling constants 

G = {g = (</ 2 , 94, 9b) G R 3 : # 2 > 0, # 4 < 0, g 6 > 0} (33) 

of the family of even sixtic potentials 

V(X,g) = g 2 X + 94X 2 + 9eX 3 , (A = z 2 ). (34) 

We refer to appendix A for the proof of the following facts: (i) the inequality 



2g 



2 >! (35) 



4 

determines an open subset of the one-cut phase G\\ (ii) the boundary Y of 
this subset is the elliptic cone 

50206 = 2gl; (36) 

(iii) the sixtic model (134"]) is singular on the curve 7 given by 

5g 2 g 6 = 2gl Ag\ = -225<? 6 2 ; (37) 

and (iv) the model is in G\ for g in V — 7 . 

We apply these results to study the deformations of ( 13"2"j) from the initial 
point g = (90, —15, 1) in V. The homogeneity in T of fl36|) implies that the 
curve g/T lies on V for all T > 0, while from ( 157)) it follows that g/T 6 7 
only at T = 60. Then g/T e Gi for all T 7^ 60 while g/60 represents the 
multicritical string model of j25[ with potential function 

VJX) = -A - ^A 2 + —A 3 . (38) 
v ; 2 4 60 K ' 

Finally, if we apply first the Bleher-Its deformation to V(A, g/T) the resulting 
coupling parameters are 

92(T,t) = (l - 1) + 4 (T,t) = |i„ go(T,t) = ^. (39) 



For the particular values g 2 = 90, g 4 = — 15, g e = 1 of (1521) a direct compu- 
tation shows that 

5 g 2 (T,t)g 6 (T,t) _ T(t - 1) 

2 g,(T,ty i+ 90 ' l4UJ 

and using (1351) we find that g(T, i) G Gi for all T 7^ 60 and i > 1. 

3. Asymptotics of the recurrence coefficients 

The main equation to determine the asymptotics of the recurrence coef- 
ficients is the discrete string equation 0] 

n 

Vz(L) n ,n-l = ^- (41) 

Here V z stands for the derivative of the potential with respect to z 

2p 

V z {z) = Y J ^g 2k z 2k ~\ (42) 
fc=i 

the subindex (n, n — 1) denotes the corresponding matrix element between 
the orthogonal polynomials 

V n (x) = P n , N ( X ), n >° ( 43 ) 

defined in (JTj), and the operator L acts on this family of polynomials as 

Lv n = v n+1 + r^jv^n-i, r 0jN = 0. (44) 
The string equation (f4T|) can be written in the form 

HA n 

^Vx(X)Un^(X) = - (A = z 2 ), (45) 



where 7 is a large positively oriented circle \X\ = R, and U n> w is the generating 
function 

U n AX) = l + 2j2(L 2k - 1 ) n ^ 1 \- k . (46) 
fc>i 



Note that L is the Lax operator of the Toda hierarchy [24j and that U n ^ is 
related to the resolvent IZ(z) = (L — 2) -1 of L by 

U n!N (X) = 1 - K{z) n , n - X - K{-z) n , n - X . (47) 

In appendix B we show that U n ^ satisfies the quadratic equation 

r n ,N {U n ,N + U n - ljN ) (U njN + U n+ltN ) = A (Ul >N - 1) . (48) 
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3.1. The continuum limit of the recurrence coefficient in the one-cut case 

We recall that our final goal is to solve the string equation ( |4"5]1 for the 
recurrence coefficient in the large N limit. If g £ G\ it has been rigorously 
established 0, E2I 27] that the asymptotics of the recurrence coefficients as 

n 

n ->■ 00, AT ->■ 00 and — ->■ T (49) 

in a neighborhood of T = 1 is given by a series 

r niJV (g)~r( e ,r,g) (e=l/JV), (50) 

of the form 



r(e,T,g) = ]Tr fc (T,g)6 2fc . (51) 

fc>0 

In particular the leading coefficient is 



a 2 



ro = T , (52) 

where (—a, a) is the eigenvalue support for the model V(X,g)/T. We write 
the asymptotics of the generating function U n ^ as a similar series 

U niN (\)~U(\,e;r), (53) 
U(X,e;r) = J2Uk^r ,...,r k )e 2k . (54) 

fc>0 

Substituting the series ( 15T1) and ( 1541) . and the corresponding shifted expan- 
sions 

r n+jiN ~ ry](e,T) = r(e,T + je), k £ Z, (55) 
C/„ +i)iV (A) ~ t%](A,e;r) = [/(A, e; r b1 ), fc £ Z, (56) 
into ( l48l . we get 



r ([/ + [/[_!]) ([/ + U {1] ) =\(U 2 -1). (57) 
Incidentally, we note as a useful consequence of ( ISTl) the linear equation 

Hi] (C/ [2] + C/[i]) ~r(U + [/ hl] ) = A ([/[!] - [/) . (58) 
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Identifying powers of e recursively in ( 15 7p or in (1581) we find that the 
coefficients Uk can be written in the form 



3A- 



where 



C7 ° = Vx^' (60) 

and the functions Uk,j(ro, . . . , r^) are polynomials of degree j in r , . . . , 
and their T derivatives. Moreover, these polynomials are homogeneous of 
degree 2k with respect to the weight w(d l T Tj) — i + 2j, and the dependence 
of Uk in rfc comes solely from 

£4,1 = 2r fc , (61) 

so that 

U k = Uo (t^T- + ■ ■ • ) > ( 62 ) 



\A-4r / 

where the dots stand for terms in r,-, and their T derivatives r', r'', . . . with 

" J J 

j = 0, . . . , k — 1. We give explicitly the polynomials £4j corresponding to 
fc = l: 



(63) 





= 2n, 




= 2r <, 


u 1>3 


= 10r (r[,) 2 



and to = 2: 

C/ 2 i = 2r 2 



^2,2 = 2^ + 2^ + 6^ + lr ri 4) , 

6 

U 2 ,s = 20ror;r; + yr r;r^ 3) + 10r 1 (r[ ) ) 2 + yr (r[ ) , ) 2 

+ 20r 1 r r[ ) , + 2r 2 rJ 4) , (64) 

U 2A = 140ro(r[ ) ) 2 r[ ) , + 56ro 2 r[ ) rS 3) + 140r 1 ro(r[ ) ) 2 + 42r 2 (r[ ) , ) 2 , 
C/ 2)5 = 924r 2 (r[ ) ) 2 ^ + 378ro(r[ ) ) 4 , 
C/ 2)6 = 2310r 2 (r[ ) ) 4 . 
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Likewise, the continuum limit of the string equation (|45p can be written 



as 

clA 



7 V A (A)[/(A,e;r)=T (65) 

Z7T1 

or in terms of the expansion coefficients U^, 

i 7T^Vx(\)U h (\; r , . . . , r k ) = 5 kfl T, k>0. (66) 

J j 27T1 

Let us introduce the function 



W(r a , g) = jf ^V A (A) = g Q w » (67) 

Then, the = equation in fl66l) is 

^(r ,g) = T, (68) 

or explicitly, 

E( 2 ;)^ 2n r« = T. (69) 

71=1 \ ' 

This is an hodograph-type equation for r$ (i.e., an equation which is linear in 
the independent variables g 2 , ■ ■ ■ , g 2p , T). Our next aim is to prove that fl66l) 
permits the recursive calculation of the r k as functions of T and g. 

For an even potential in the one-cut case, using the variable A = z 2 and 
the relation between a and r given in equation (|51Z|) . we find that Wi(z) = 
V z 2 — a 2 = a/A — Atq. Therefore the definition (1171) of the polynomial h(X) 
can be written as 

2 V\V X (X) =h ( X \ + t x -i^ a^oo. (70) 
V A - 4r 

Consequently 

°W(r , e ) = 2<fl± /f> =h(X) . (71) 

or / 7 27ri (A - Ar y/ 2 A=4r 

Thus, given g G G\ equation ( [71]) implies that d ro W(r , g ) 7^ 0, and the 
hodograph equation ( |68i) defines implicitly tq as a locally smooth function of 
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T and g in a neighborhood of T = 1 and of go- The remaining equations (I66j) 
can be written as 



3A: 

^ / W j (r ,g)U kJ (r ,...,r k ) = 0, k>l, (72) 

where 

^(r ,g) - £ ^i (A _ 4ro)i+ : " Wj -I)!!" (73) 

From flBTj) and f!73|) it follows that the term (<9 ro Jy(r , g))r^ in equations f J72|) 
equals a sum of terms of in r , . . . , r k -\ and their T derivatives, and therefore 
define recursively the coefficients as locally smooth functions of T and g 
in a neighborhood of To = 1 and of go. 

We recall again that the coefficients [4j(r , . . . , r>) in (172"]) are polynomi- 
als of degree j in Tq, . . . , r& and their T derivatives. Repeated differentiation 
of the hodograph equation (16H|) with respect to T give the T derivatives of 
r as a rational function of r : 

1 n d? o W 

r ' = d^w' r ° = ~(^p ••• (74) 

and we can effectively solve equations (1721) for as a rational function of tq. 
Using the standard notation W', W", . . . , W^> for the derivatives of W with 
respect to r we find 

_ 2 (in 2 - 

ri - r ° i2(wr ' 

and 

X + tqY 

r2 = r °ii40(^r' (76) 



where 



X = 700W' (W"f - 910 (W'f (W"f W'" + 118 (W'f (W"'f 

+ 180 (W'f W"W {4) - 18 (W) 4 W {5 \ (77) 

F = -980(^ ,/ ) 5 + 1760^ , (W r// ) 3 W AW - 545(W r/ ) 2 ^ // (^ w ) 2 
- 420 (W'f (W'f W {A) + 102 (W'f W"'W (4) 
+ QA(W'fW"W^-h(W'fW^\ (78) 
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which in turn show that the are rational functions of r . These expressions 
for r\ and agree with Eq. (4.25) and Eq. (4.56) obtained by a different 



method in 19 



Finally, we remark that our method to calculate the coefficients of the 
large N expansion ultimately depends only on equations ( 15"7|) and (JB"5|) . which 
are invariant under the symmetry transformation 

(e~T,g) = -(e,T,g), c > 0. (79) 
c 

Hence it follows that r(e, T, g) = r(e, T, g) and consequently 

r*(r,g) = ^r fc (l,g/r), fc>0. (80) 

TTie recurrence coefficient under the Bleher-Its deformation in the reg- 
ular one-cut case 

Let g G Gi such that its Bleher-Its deformation is in the regular one-cut 
case, i.e., g(£) G G\ for all £ > 1. Then we can apply the results of the 
previous subsection with g replaced by g(£) to conclude that in the limit (jHJ) 

r n)JV (g(£))~r(e,T,g(£)) (81) 

where the coefficients r k (T, g(£)) of the asymptotic series 

r(6,T,g(£)) = 5> fc (T,g(£))e 2fc (82) 

are determined by 

W(r ,g(t)) = T, (83) 

3fe 

^H/,(r ,g(£))f/ fcj (r ,...,r fc ) = 0, fc > 1, (84) 

3=1 

as smooth functions of t and T for £ > 1 and T near T = 1. Our aim 
in this subsection is to find a reformulation of fl83l) - fl84l) that decouples the 
dependence on t and g. 

The string equation ( 1651) for the deformed model is 

i4y A (A,g(£))[/(A,e;r)=r. (85) 

Z7T1 
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If we substitute in this equation 

V x (\,g(t)) = l-~ + ~V x (j, g y (86) 
and take into account that U ~l + 2r/AasA — > oo, we find that 

2(t - 1) 7 + i^(T' g ) C,(A ' f,r) = T ' (87) 

or with the change of variable A — > Xt 

2(*-l)7+ <f^:V x (\,g)U(\t,e,r)=T. (88) 

t J y ZTil 

Note that the generating function U(X, e;r) is uniquely determined by (|57l) 
and the asymptotic behavior U(X,e;r) ~ 1 as A — > oo. Since U(Xt, e; r) 
satisfies ( 1571) with the substitution r — > r/t and U(\t,e;r) ~ 1 as A — > oo, 
we conclude that £/(At, e; r) = £/(A, e; r/t). Alternatively, we can arrive at 
the same conclusion directly from the explicit expressions ( 159]) and fl60|) for 
the Uk and from the fact that the functions Ukj are polynomials of degree j 
in 7*0, ... , Tfc and their T derivatives. Therefore, if we denote 

t(e, T, t, g) = r(£ ' T ; g(t)) = £ *> s)* 5 *, ( 89 ) 

fc>0 

equation flHHl) becomes 

2(t - l)r + $ ^^a(A, g)C/(A, e, r) = T, (90) 

or equivalently 

2(t-l)t + W(t ,g) = T, (91) 

2(t-l)t fe + ^^(to,g)£4j(t ,...,r fe ) = 0, fc>l. (92) 
i=i 

Note that the only changes introduced by the Bleher-Its deformation in our 
calculation of the recurrence coefficient are the first term 2(t — l)t& and the 
substitution r*i — > tj in Ukj- 
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4. Topological expansions in the regular one-cut case 

In this section we implement our method to calculate the coefficients 
of the topological expansion by means of our equations (l9~Tj) - (l92l and the 
Bleher-Its representation of the free energy, which we repeat here for conve- 
nience: 



F N (g) = F$ + 



1 - t 



r N>N (g{t)) (rjv_i,Ar(g(t)) + r N+1>N {g(t))) - - 



dt. 



(93) 

In the first subsection we derive integral expressions for the coefficients of 
the topological expansion. Next, drawing on an idea that Bleher and Its 0] 
used to calculate the leading term F(°)(g) of the expansion, we give an ef- 
ficient procedure to calculate higher-order coefficients and compute explicit 
expressions of the first four coefficients for general models. Finally, we apply 
our results to three widely studied models. 

4-1. Expressions for the coefficients of the topological expansion 

Equation (I9"3l involves only the recurrence coefficients r N>N and r N ± 1>N . 
Therefore to study the limit (jSj) in the regular one-cut case we need only the 
asymptotic series for r n ^ in a neighborhood of T = 1: 

r N As(t)) ~ r(e, 1, g(t)) = tt(e, 1, t, g), (94) 

r N ±i,N ~ T{e, 1 ± e, g(t)) = tt(e, 1 ± e, t, g). (95) 
Substituting these expansions in the Bleher-Its formula we have 

/oo 
(l-t)/(e,t,g)dt, (96) 

where (using again our shifting notation t[±i](e, 1, t, g) = r(e, 1 ± e, t, g)) 

/(e,t,g) =t(e,l,t,g) (r hl] (e,l,t,g) + r [1] (e,l,t,g)) - ^. (97) 

With the method discussed in the previous section we can readily obtain an 
expansion 

/M,g) = ]T/ fe (i,g)e 2fc , (98) 

fc>0 
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where the first five coefficients are 




(99) 
(100) 
(101) 




(102) 



(103) 



Here the primes denote derivatives with respect to T evaluated at T = 1. 
Note also that the term /& has weight w(fk) = 2k. The analysis of Bleher and 
Its in [2] shows that in the regular one-cut case it is legitimate to perform term 
by term integration in ( |96|) . which yields the following topological expansion 
of the free energy: 



Therefore the direct method to calculate the coefficients of the topological 
expansion (11041) is as follows: first, use equations (|9TT) -(l92l) to determine the 
coefficients tk, then use equation ( |9"7|) to find the and finally perform the 
integration with respect to t in equation fllOSp . 

4-2. Efficient calculation of the coefficients of the topological expansion 

The direct method to determine the coefficients of the topological expan- 
sion outlined in the preceding paragraph requires explicit calculation of the 
functions tk(T, t, g) which is, except in the simplest difficult task. 




(104) 



where 




(105) 
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In [2| Bleher and Its used an ingenious idea to determine the leading coef- 
ficient F^(g) for general models V(\,g). In this section we will show that 
the same idea can be applied to evaluate higher order coefficients F^'(g). 

It follows from our previous results that the functions fk{t,g) can be 
written as rational functions of t and to(t, g). Now, if we denote 



£ = tb(l,*,g), 
the hodograph equation fl9Tj) at T = 1 implies that 

t = l + i(l-W(fcg)), 



(106) 



(107) 



which suggests to use £ as integration variable in (I105p . At the lower limit of 
integration t = 1 we have that £ = rn(l, g), while at the upper limit the new 
variable £ ~ 1 / (2t) — >■ as t — > 00 [2] . Therefore, with a trivial sign change 
absorbed in the definition of Rk, equation (I105P can be written as 



2r(*)| 



g) 



i? fc (e, g )de, fc>i, 



(108) 



where the -Rfc(£, g) are rational functions of £. Note that the evaluation of 
these integrals yields the coefficients as functions of g and r = r (l, g). 

The presence of the term l/(2£ 2 ) in the expression of fo requires a slightly 
different integration process [2j to calculate F(°)(g). Using equations (|99|) . 
(DUD and (JTU7D w e find: 



lim 

T— ^OO 



lim 

t— Kxi 



1/2t 



;i-t(0)e 2 dt(0 + olnr-- 



l/2r 1 

+ -lnr 

ro 2 



+ 



''() 







1 

+ 2 



1 

" 2 

(W(0 - 1) W'(fld£. 



(109) 



Hence, taking into account that W(ro, 

In ro 



1 at T = 1, it follows that 



F (0) (g) 



ln2 3 

~-4 + 



''() 



w — w 2 

2 



110) 
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Note that (since W does not have a constant term) the integrand in the last 
term of fill 01) is a polynomial in £. 

Let us proceed now with the calculation of F^(g). We repeat here for 
convenience equation (l9Tj) and particularize equation fl92|) for k — 1: 

2(t-l)tb + W(tb,g)=T, 



2(t - l)ti + ^ W^(tb, g)t/i,,(t , ti) = 0. 

Differentiating with respect to T in fillip it follows that 

1 „ 3W, 



111) 

;ii2) 



u 2(f-l + Wi)' u 2(t- 1 + W^i) 3 : 
and substituting these expressions in (j!12p we find 



(113) 



^i — 
Therefore 

fx = t (4ti + to) = t 



3W% 5W 3 



2{t - 1 + W^) 4 4(t - 1 + Wi 



_ 10W 3 t + 3W 2 
(t-l + Wx) A ~ 2(t - 1 + W^ 3 



(114) 



(115) 



and changing the variable from t to £ with ( 11071) we obtain the following 
expression for fW(g): 



F«(g) 



(W-l) 



£ 2 24Wf 



10^3 + 3iy 2 



;i + 2fWi - Wf (1 + 2^i - W) 2 _ 
Using f J73|) we can identify the integrand as a total derivative: 
f (-1 + W)W (-3 + - 3fW + 2£ 2 W") £ 2 (-l + W)W 



£d£. (116) 



12 d£ 

Thus we arrive at the result 



12(1 -W + ZW) 6 

( In {-£W' + W - I) 



12(1- W + iW) 
(W - l)j 2 W" 

'-ZW + w-i) 



F«(g) = ^ln(r Vy'(ro,g)) 



(117) 



118) 
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in agreement with Eq. (4.27) of 20 . 

The calculation of the next coefficient (g) of the topological expansion 
is entirely similar, and we omit the intermediate steps which are easily per- 
formed with a symbolic computation program. Using (191 p . fl92l) for k — 1,2 
and (llOip we find 



R 2 (C,g) = 3(1 -W) 



■56448W,' 



f 7 

- 2j (103681U 2 4 + 8M80£W%W 3 ) 
+ [420W 2 3 + 10800^ 2 2 W 3 



+ £ 2 (21800WW 3 2 + 23520 W%W 4 )] 

^ [260W 2 W 3 + £ (lllOWf + 2380W 2 Wa) 



+ £ 2 (476OW3W4 + 5376W 2 W £ 

+ ^ (35W 4 + 336f W 5 + 770£ 2 W 6 ) 
cr 4 



where 



a = l-W + 2£Wi = 1-W + £W'. 
Using (175)) we can again identify i2 2 (6j g) as a total derivative, 



R2& 



2880 d£ 



a' 



■280(-l + 



tf\4 



£ 6 



+ -^(-1 + W) (300(W) 3 + 400£(W") 2 W'") 

:5 



1 X 



-^[5Qaw"f + (-i + wy. 

(260W"W" + 58£(W'") 2 + 88£W"W {4) )] 

+ [-9(W") 2 + b8£W"W"' 

+ (-1 + W) (361U (4) + 10£IU {5) )] 

+ \ (-12W" + A£W"' - lO^W^) - 
a 6 o l 



(119) 
(120) 



;i2i) 
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and we finally obtain 



+ 



+ 



7r W"(r f 



+ 



+ 



240 240r 2 W'(r ) 2 360W(r ) 5 
W"(r ) (9W"(r ) - 58r W^(r )) 

2880^'(r ) 4 
6W"(r ) - 2r W^(r ) + br 2 W^{r Q ) 
1440r W // (r ) 3 



(122) 



This expression agrees with the result of [20] (up to a trivial mistake in the 
sign of the first fraction of his Eq. (1.18), which is fixed in his application to 
the quartic model in Eq. (4.59)). Our method can be easily carried on further 
with a symbolic computation program. For example, the next coefficient 
turns out to be 



1 



W"{r ) 



1008 



1008r 4 l¥'(r ) 4 
' [ 15W"{r ) 2 - AW {3) {r )W'{r ) 



6048r 2 W'{r y 
1 



15W"(r ) 3 + W w {r )W'{r )' 



Q048r W'(r y 

- 10W {3 \r )W(r )W"(r ) 

1 

725760W(r ) 8 . 
+ 2001^ (3) (r ) 2 W /, (r ) 2 + 300W (4) (r )W'(r ) 2 W"(r ) 

- 1800W {3) (r Q )W'(ro)W"(r )^ 



1575W'(r ) 4 - 24W {5) (r )W'(r y 



r 



362880M/'(r ) 9 
+ lQUW {b) {r )W'{r G ) 3 W"{r^ 
+ 2A88W i3) (r )W w (r )W'(r ) 3 

- 10170W r(4 Vo)W r '(ro) 2 W r// (r ) 2 
+ 40110iy (3) (ro)W /, (r )W /,/ (ro) 3 

- 12783W {3) (r ) 2 W'(r ) 2 W"(r ) 



- 21420 W"{r Q f - 133W^(r )W'(r o y 



362880W(r ) 



10 



34300Vr"(ro) 6 - 35W (7) (r )W(r ) £ 
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+ 607W w (r ) 2 W'(r y - 2915^ {3) (r ) 3 H // (r ) 3 
+ 539W {6) (r )W'(r ) 4 W"(r ) 
+ 1006W {3) (r )W {5) (r )W f (r ) 4 

- A28AW^(r )W'(r ) 3 W"{r ) 2 
+ 22260W i4) (r )W'(r ) 2 W"(r ) 3 

- %l$mW {3 \r )W'{r )W"{r Q y 
+ 43050W {3) (r ) 2 W'(r ) 2 W"(r ) 2 

- lUb2W {3 \r )W {4 \r Q )W\r Q ) 3 W'\r )\. (123) 

This expression reduces to the result found by Shirokura except for the re- 
placement of 388 by 300 in the fourth coefficient of Eq in Eq. (50) of 19 



We apply now flllOp . (11181) and (11221) to obtain explicit expressions for 
the first three coefficients of the quartic, two-valence and sixtic models (we 
omit the lengthy expressions for F^(g) which are obtained in exactly the 
same way using (11231) ). 

4-2.1. The quartic model in the regular one- cut case 

We have shown in section [2~T1 that if (#2, #4) G Gi the quartic model (1231 
is in the regular one-cut case. Using equations (jllOp . f II 1 8 j) and fll 22j) with 

W(r , g) = 2g 2 r + 12g 4 r 2 (124) 

we find 

^ (0) (g) = -§ + I (92r ) - I (92r ) 2 - ^ In (2r ) , (125) 

F( 1 '(g) = lln(2(l- W )) ) (126) 

p(2 )/„n (2^ 2 ro-l) 3 (41 + 21^r -6(^r ) 2 ) 

F (g) = H520(l -g 2 r Q f ' (12?) 

where r is the positive root of the hodograph equation (|69|) 

2g 2 r + 12 9i r 2 = 1, (128) 

namely 



r = ~ g2 + ^ + 12g ' . (129) 
Lzg 4 
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4-2.2. Two-valence models 
For the two-valence models 

V(\%) = g 2 \ + g 2u \\ v>2, (130) 

in the region g 2 > 0, g 2u > 0, the Bleher-Its deformed potential is a convex 
function of z for all t > 1, and therefore these models are in the regular 
one-cut case. Using 

W(r ,g) = 2g 2 r + u^\g 2u r v (131) 

we find: 

^. N 3(i/-l) (z/-l)(2i/ + l), N (^-l) 2 , „ 1, , N 
4z/ i/(i/ + 1) v[y + 1) 2 

(132) 

F( 1 )(g) = lln(z/-(z/-l)2^ 2 r ), (133) 

F w (e) = gwJH x 

vs; 2880(i/ - 2(i/ - l)# 2 r ) 5 

- z/ 3 (8z/ 2 + 5z/ - 1) + 2v 2 {y - l)(16z/ 2 + 40z/ - l)(g 2 r ) 

- Av{v - l) 2 (8z/ 2 -v + 44)(# 2 r ) 2 

- 96(i/ - l) 3 (4z/ + l)(# 2 r ) 3 + 192(i/ - l) 4 (^ 2 r ) 4 ] • (134) 
Here r is determined as the positive root of the hodograph equation 

2^0 + 1/^^ = 1. (135) 

Note that f ll32l) - fll34p reduce to the quartic results for v — 2. For 02 = 1/2 



the expressions (I132l) - fll34p reduce to the corresponding results in [22|, [28 



4-2.3. Sixtic models 

Our last application concerns sixtic potentials 



V{\^) = g 2 \ + g A \ 2 + g^ (136) 
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in the region g 2 > 0, g 4 > 0, g 6 > 0. Again the convexity argument proves 
that these models are in the regular one-cut case. Using 



we find: 



W(r , g) = 2g 2 r + 12g,r 2 + 60g 6 r 3 (137) 



^ (0) (g) = " g + 1 (92r ) - 3 (g 2 r y + 1 fo 4 r 2 ) - jj (<? 4 r 2 ) 2 - || Wq 3 - - In (2r ) , 

(138) 

F«(g) = ^ In (3 - % 2 r - 12^) , (139) 

^f>\ /\ 1 593 

^ () (g) = -7T77T + 



+ 



240 720(122 4 r 2 + 43 2 ro-3y 
169# 2 + 2928^4 - 1716^ 2 r 



+ 



+ 



720^ (l2g A r 2 + 4g 2 r - 3) 3 

224# 2 4 + 7587g4g| + 45765g| - (57888<& 2 + 6756^1% 

6480g 2 (12g 4 r 2 + 4g 2 r -3) 4 
7 (6gj + 81^1 + 243g| - (8g 5 2 + 126^1 + 486^ 2 )r ) 
405g 2 (12g 4 r 2 + 4g 2 r -3f 

(140) 

In this case Tq is the positive root of the hodograph equation 

2(7 2 r + 12(7 4 r 2 + 60^r 3 = l. (141) 



5. Counting numbers 

The methods of Ercolani and McLaughlin Q prove the existence of the 
topological expansion for matrix models 

WA) = ^ + X>^, (142) 

3=1 

under the hypothesis that there exists a path in the space of coupling con- 
stants connecting t = (t 2 ,t 4 , . . . ,t 2u ) to the origin 0. The corresponding 
coefficients F^ k \t) are analytic functions of t near the origin and their Tay- 
lor expansions determine graphical enumeration numbers. 
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To put these results in context, we briefly recall that a /c-map is a graph 
which is embedded into a surface of genus k in such a way that (i) the edges 
do not intersect and (ii) dissecting the surface along the edges decomposes 
it into faces which arc homcomorphic to a disk. We can formulate the result 
of j3] as the following representation: if 

j?W( t ) = - -r^ r(-t 2 ) n2 • • • (-t 2p r^ k (n 2 , . . .,n 2p ), (143) 

then /tfc(^2, • • • , n 2 p) is the number of connected fc-maps with a number n 2 j of 
2j-valent vertices in which all the vertices are labeled as distinct and all the 
edges emanating from each vertex are labelled as distinct as well 

It is straightforward to rephrase our results of the previous section in 
the notation of (11431) and therefore we have a direct method to calculate 
the Kk(n 2 , . . . ,n 2p ): we introduce the change of variable A' = A/2 to reduce 
Vbm(A) to the form (1201) . which in turn implies the following relation between 
our set of g coupling constants and t: 

g 2k (t) = 5 K2 + 2 k t 2k} k = l,...,p. (144) 

As an application of this first, direct procedure, we have carried out these 
substitutions in the topological expansion of the quartic model (j!25p - (|127|) . 
expanded in Taylor series these coefficients, and found the corresponding 
counting numbers Kk(n 2 ,n4), the first of which we present in table [TJ 

However, since the counting numbers depend on the coefficients of the 
Taylor expansion of the F^ k \ drawing on our analysis of the previous sec- 
tion we can calculate directly these Taylor coefficients without requiring the 
explicit evaluation of the F^ k ' themselves, which turns out to be much more 
efficient. In essence, the idea is to obtain first the Taylor expansion in t of 
the integrand of f 1 1 8 j) . and subsequently to perform the integration in (1105)) 
with respect to the Bleher-Its deformation parameter t term by term. A sim- 
ilar approach has been recently used in ji| to investigate counting numbers 
in the cubic model. In detail, the procedure is as follows: 

1. Use ( 1T44D to write Eq. (ED at T = 1 as 

j=l N ' 
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Table 1: Lowest three counting numbers Kk(ri2, n±) for the quartic model V(X) = 
and all combinations of the up to 4. 
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and then use implicit differentiation to calculate the Taylor expansion 
at t = 

3i,---,3 P > 

2. Use the string equations fl§T|) - fl92"j) to write fk as rational functions of 
to- Then substitute (11461) in the resulting expressions and determine 
the Taylor expansion of fk at t = 0. 

3. Perform the integration 

/oo 
(l-t)f k (t,t)dt, (147) 

term by term and find the numbers Hk{n 2l ■ ■ ■ , n 2p)- 

We have implemented this strategy for the two-valence and for the sixtic 
models, and present some of our results in table [2j The /%(ri2, . . . , n 2v ) 
grow quickly in number and in magnitude, and we give a complete table 
only up to n 8 = 2. Note that the results in table [2] with = agree 
with the corresponding results in table [TJ Our results also agree with those 



for Ki(n 2} 0, . . . , n 2u ) and K 2 (n 2} 0, . . . , n 2u ) in |28j. Similarly, in table [3] we 
present all the fifth nonvanishing counting numbers K4(n 2 , 714, Uq) with < 
rii < 3 for the sixtic model. 



6. Singular one-cut cases 

Let g G G be such that its Bleher-Its deformation is in the singular one- 
cut case, i.e., g(t) G G\ for all t > 1 but g = g(l) determines a singular 
model. Our discussion in sections [3] and H] shows that the integrals for the 
coefficient of the topological expansion (jl04p — (1105p converge provided that 
the function to(T, t, g) is smooth near (T, t) = (1, 1). For this type of singular 
cases the topological expansion f llU4j) (I105j) still exists. An example of this 
situation is the quartic model ([231 for g 2 = —2^fgl. 

However, in the singular one-cut cases where the function to(T, t, g) is 
not smooth near (T,t) = (1,1) the topological expansion (I104p - (I105I) is ill- 
defined, because the integrals defining the coefficients for k > 1 diverge. 
This is the case of the Brezin-Marinari-Parisi model fl38l) . We next discuss 
this critical behavior in general and present a method of regularization. 
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Tabic 2: Lowest four counting numbers Kk(n2 1 ni 1 no) for the sixtic model V(X) = 52A + 
54A 2 + geX 3 and all combinations of the rn up to 2. 
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Tabic 3: Nonvanishing counting numbers Ki(n2, rn, ne) with < rii < 3 for the sixtic 
model V(X) = g 2 X + .94A 2 + g 6 X 3 . 



n 2 


n 4 


n 6 


K 4 (n 2 ,n 4 ,n 6 ) 





1 


3 


1143525600 





2 


3 


2201217638400 





3 


2 


24069830400 





3 


3 


2836746385920000 


1 


1 


3 


25157563200 


1 


2 


3 


57231658598400 


1 


3 


2 


577675929600 


1 


3 


3 


85102391577600000 


2 


1 


3 


603781516800 


2 


2 


3 


1602486440755200 


2 


3 


2 


15019574169600 


2 


3 


3 


2723276530483200000 


3 


1 


3 


15698319436800 


3 


2 


3 


48074593222656000 


3 


3 


2 


420548076748800 


3 


3 


3 


92591402036428800000 
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6.1. Critical behavior and a triple- scaling method of regularization 

Let us consider a singular one-cut deformation such that the string equa- 
tion (l65l) has a critical point of order m > 2 at (tq,T) = (r c , 1). That it to 

say, 

W(r c ,g) = l, (148) 

d ro W(r c , g) = • ■ ■ = d™" 1 W(r c , g) = 0, (r c , g) + 0. (149) 

In this case the implicit function theorem does not apply to equation ( 19 ip 
near (T , t ) = (1, 1) with to(l, 1, g) = r c . In fact, Xq(T, 1, g) can be expanded 
m powers of (T - lf/ m and, provided that r c + 0, r (l,t,g) can be also 
expanded in powers of (t - l) 1/m . As a consequence, the system f l9~Tj) - f[9"2"j) 
does not yield an appropriate asymptotic series to generate the free-energy 
expansion (|96|) . In order to regularize this critical behavior it is natural to 
introduce two scaling variables x and y in the form 

T = 1 + e m x, t = l + e m y, (150) 

where 



1 \ 2m + l 

e^+i = [ — ) . (151) 



In terms of these scaled variables the string equation (|90|) reads 

2e m yr + j> |Va(A, g)C/(A, e; t) = 1 + e™x, (152) 
and we will prove now that there are solutions of the form 

t(e, x, y,g) = r c + J2 * [k] ( x i Vi g)e fc - (153) 

fe>i 

Note that the shifts T — > T ± e correspond to x — > x ± e 1 / 2 , and therefore 
Z7(A, e; t) is determined by the quadratic equation 

r (U + Upq) (U + U m ) = A [U 2 - l) , (154) 

where we have denoted /rw(^) = /(a? + fce 1 / 2 ). 

The corresponding expansion of U to be substituted in the string equa- 
tion fTl52|) is 

C/(A,e) = ^f/ [fc] (A;r [1] ,...,r [fc] )e fc (155) 

fc>0 
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where 



c 



From (1154ft it follows that the U^'^ are polynomials in tj^K . . . , r' fc_ - ?+ ^ and 
their x derivatives, which can be determined recursively |30J]. In particular 

U [k ' l] = 2x [k] . (158) 
The first few of these coefficients are 

u{ 2,2] = 6(^1)2 + 2^ (159) 

^ = 12^^ + 2^^1+2^1 + 1^, (160) 

^[3,3] = 20(rW) 3 + 10r c (rW) 2 + 20r c rWrW+2r c 2 4 1 ] a;x . (161) 
As a consequence of the quadratic equation we find the linear equation 

t[i] (U m + U m ) -x(U + U FT] )=\ (U m - U) , (162) 
which in turn leads immediately to the recursion relation 

dxU [k+iMi] = ( rc(9 3 + 4x [i] dx + 2r [l]) u^ k \ C/[°-°] = 1. (163) 

This relation implies that the U^ k,k \x^) are the well-known Gel'fand-Dikii 
differential polynomials of the KdV theory (3l| . 

We now substitute f ll 55j) into (11521) and take into account (1731) . (I148j) 
and (1X491 to obtain 

2e m yt + £ W^(r c , g)tf W (e, r) = ft. (164) 

fc>m 

Finally, collecting powers of e we find the system of equations: 

W m (r c ,g)U^ m \^)=x-2r c y, (165) 

m+k 

2yx [k] + Wj(r c , g)U^' m+k] (x [1] , t [m+k ~ j+1] ) = 0, k > 1. (166) 

32 



The first equation constrains y, g) to be of the form u(x — 2r c y), with 

u(x) being a solution of the m-th member of the Painleve I hierarchy. The 
subsequent equations (I166P give for each coefficient xP^(x,y) with (k > 2) an 
ordinary differential equation in the x variable involving the previous coeffi- 
cients (1 < j < k). The characterization of the appropriate solutions of 
these ordinary differential equations is a difficult problem deeply connected 
to the regularization of the free energy expansion jij. 

To regularize the expression ( 196]) we partition the domain [1, +oo) of the t 
variable into an inner region [1, l+<5(e)] and an outer region [l+#(e), +oo). In 
the inner region we assume the triple-scaling limit asymptotics fl 1 5 3 [) for the 
recurrence coefficient, while we assume the regular one-cut asymptotics 
in the outer region. Thus, we have 



F N (g)-F^I 1 (e,g) + I 2 (e,g), N^oo, (167) 



where 



h(e,g) = -e 2m f m/C yf 1 \e,y,g)dy, (168) 
Jo 



f (1 \t,y,g) = t(e,0,y,g) (tpi](e,0,y, g) + t^e, 0, y, g)) 
and 



2(l + e m y) 2 ' 

(169) 



J 2 (e,g)= / (l-t)/( 2 )(e,l,t,g)dt, (170) 

/ (2) (e,t,g) =t(e,l,t,g) (t[_i](e, 1, t, g) +t {l] (e, 1, t, g)) - — . (171) 

To ensure that the result is independent of the choice of 5(e), the asymptotic 
series /^(e, y,g) and f( 2 \e,t,g) must be matched on some appropriate in- 
termediate region overlapping the inner and outer regions. It is at this point 
where the conditions to determine the coefficients of (11531) emerge. 

6.2. The Brezin-Marinari-Parisi critical model 

We illustrate these ideas with the Brezin-Marinari-Parisi critical poten- 
tial (1381) . which belongs to the singular one-cut case under the Bleher-Its 
deformation. In this case 

W(x 0l g)=x 3 -3x 2 + 3x 0} (172) 
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where g = (g^g^ge) = (3/2,-1/4,1/60), and according to ffT4^D- ffT4^j) we 
have a critical point of order m = 3 at r c = 1. Thus, from (I165p it follows 
that = u(x — 2y), where u(x) is a solution of the second member of the 
Painleve I hierarchy 

Uxxxx + 10uu xx + 5u x + 10m 3 = lOx. (173) 

To perform the matching between the triple-scaling and the one-cut reg- 
ular asymptotics, we note that as t — > 1 + , from fl9T|) and (1921) we have 

r (l,t,g)~l-2 1 / 3 (t-l)i, (174) 

r' '(l,t,g)~3- 2 2- 2 / 3 (t-l)-f, (175) 

Hence we get 

f (2) (e, t, g) ~ ( 2 r 2 e\(^ + tg) ~ \ - a^V* (177) 

Likewise, in the inner region we have 

/ (1) (e,y,g)~| + e4tl 1 ](0,y,g). (178) 

In the matching region we must have both t — > 1 + and y — >■ +oo as N — > oo. 
Therefore, the matching between (I177P and (11781) is achieved provided u(x) 
is a solution of (I173p such that 

u(-2y) ~ -2 1/3 y 1/3 , ?/ -»■ +oo. (179) 

This asymptotic behavior determines a unique formal expansion of the form 
x 1 / 3 times a series in powers of x~ 7//3 , which solves (11731) . 



7. Concluding remarks 

In this paper we have developed a method to compute the large N ex- 
pansion of the free energy of Hermitian matrix models (jlj) from the large TV 
expansion of the recurrence coefficients of the associated family of orthogonal 
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polynomials. It is based on the Bleher-Its deformation, on its associated inte- 
gral representation of the free energy, and on a method for solving the string 
equation which uses the resolvent of the Lax operator of the underlying Toda 
hierarchy. Combining these ingredients we provide a procedure, suitable for 
symbolic computation, to characterize the structure of the coefficients F^(g) 
of the topological expansion of the free energy. The procedure can be also 
used efficiently for the explicit evaluation of these coefficients. As an illus- 
trative application we compute the expressions of f(*'(g) (ife = 0,...,3) for 
general matrix models and check their agreement with the expressions derived 
using the Euler-Maclaurin summation formula in the Bessis-Itzykson-Zuber 
method. 

The main application of our study is a convenient method to compute 
generating functions for the enumeration of labeled fc-maps. It relies on 
the structure of the integral representations of the coefficients F^ k \g) and 
does not require the explicit expressions of these coefficients. We apply this 
method to elaborate several tables of numbers of fc-maps with two and three 
valences and up to genus k — 4. 

Finally, in order to illustrate the regularization of singular models within 
our scheme we have formulated a triple-scaling method to regularize singular 
one-cut models. 

Although in this paper we restricted our analysis to the genus expansions 
of one-cut even models, since both the Bleher-Its representation of the free- 
energy [2| and the method for solving the string equation using the resolvent 
of the Lax operator 3(]| can be used for models associated with general (not 
necessarily even) potentials, there is no obstacle to apply our analysis to 
these problems too. Furthermore, Bleher and Its [2] applied the integral rep- 
resentation (fTUj) to determine a three-term large N asymptotic expansion for 
the free energy for the quartic model in the neighborhood of a critical point 
at the boundary between the phases G± and Gi- The third of these terms, 



which involves the Tracy- Widom distribution function [32], represents a non 



perturbative effect. Thus, it is plausible to generalize the present scheme to 
characterize nonperturbative contributions to the large N asymptotics of the 
free energy of multi-cut models. This generalization would require the use of 
nonperturbative solutions of the string equation, e.g., as in the trans-series 



method considered by Marino 33]. In particular it would be interesting to 



investigate the asymptotic behavior of the free energy in critical processes 



such as the birth of a cut in the eigenvalue support [34j, |35|, |36 
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Appendix A: Multi-cut models 



We first recall the following upper bound [36|] for the number of cuts q of 
a model with the potential V(z): 

degV(z) 

q<p= 2 • (180) 

The conditions that determine the actual value of q among those allowed by 
this bound can be stated in terms of the function h(z) defined in ( ITT)) . It can 
be shown [lo[ that in the g-cut case: 

1. The endpoints of J satisfy the equations 

"Otj+l 



h(x)w lj+ (x)dx = 0, j = l,...,q- 1, (181) 

^;MicU = 0, j = 0,...,g-l, (182) 

where 7 is a large positively oriented loop around J. Moreover, since 
fj p(x)dx = 1 we must have 

I h{z)w x (z)dz = -Am. (183) 

J-y 

2. The following inequalities hold: 

/ h^w^dx' < 0, for x < a u (184) 

J X 

f x 

/ h{x')w\{x')dx > 0, for [3j < x < aj+i, j = 1, . . . q — 1, (185) 
J Pi 

[ /i(i>i(x')ds' > 0, foix>P q . (186) 

J0c 
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Equations (I18ip - (ll83j) are 2q conditions that the 2q unknowns ai,...,(3 q 
must satisfy. However, these equations may not be sufficient to determine 
uniquely q because they may have admissible solutions for different values of 
q. If this is the case, the additional condition p(x) > for all x G J and the 
inequalities fll84[) — fll86j) characterize uniquely the solution of the problem. A 
model is said to be a regular if h(x) ^ on J and the inequalities ( I184p -( 1IS5|) 
are strict. Otherwise it is called singular. 

Sixtic potentials 

Let us consider the family of sixtic potentials 

V{X) = g 2 \ + g,\ 2 + g & \\ (A = z 2 ), (187) 

in the region of coupling constants 

G = {g = {g 2 , <? 4 , <7e) e R 3 : 92 > 0, g, < 0, g 6 > 0}. (188) 

For q = 1 and J = (—a, a) we have 

{— \x 2 — a 2 ! 1 / 2 for £ < —a, 
i|x 2 — a 2 ] 1 / 2 for —a < x < a, (189) 
|x 2 - a 2 1 1/2 for x > a, 

and 

9 

/i(x) = 6sf 6 x 4 + (45(4 + 3sf 6 a 2 )a; 2 + -g 6 a 4 + 2g^a 2 + 2# 2 . (190) 



Equation (jzzj) reads 

15# 6 A 3 + 12# 4 ^ 2 + - 16 = 0, (A = a 2 ). (191) 
Completing squares in the expression of /i(x) we have 



*m = «*^ + t + £J + +2 -"C- (192) 

Hence the function h(x) is strictly positive for all x G R provided that 

> I- (193) 

2 #4 
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Thus, the inequalities (1184p - (ll86p are strictly satisfied and p(x) > on J. 
Moreover, since the critical points of the polynomial in the left-hand side 
of (pUD are 

then (j!93p implies that there exists a unique positive solution A of (j!9ip . 
Therefore (11931) determines an open subset of G% with boundary T given by 

2gj = 5g 2 g 6 . (195) 

Given g G T it follows from HI 92[) that the function /i(x) is strictly positive 
on R unless 

c? = (196) 
155-6 

in which h(x) vanishes at x — ±a. But this value of a 2 satisfies equa- 
tion (Jill]) only if 

4^ 4 3 = -225<? 6 2 . (197) 
Hence, along the curve 7 given by 

2gj = hg 2 g & , Ag\ = -225^, (198) 

the model is singular because h(x) vanishes at the end-points ±a of the 
eigenvalue support, whereas for points in V — 7 the model is in G\. 

Appendix B: The quadratic equation for U n ^ 

For clarity in this appendix we drop the subindex iV from U n ^. Thus, 
consider the function 

U n (X) = l + 2 J2(L 2k - l ) n ,n-i\- k i (199) 

fc>l 

where the matrix elements (L x ) n>n _i are calculated in the basis of orthog- 
onal polynomials 

v n (x) = P n , N {x), n > 0. (200) 
From ([6]) it follows that L 2k ~ l v n (x) = x 2k ~ 1 v n (x), and therefore 

-2k-l\ 1 / „2k-l„ 



(Zr )„ >n _i = / x ife - 1 t;„(x)t; n _ 1 (x)d/i(x), (201) 
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where 

d/i(x) = e~ NV{x) dx. (202) 

Hence, 

_2 f°° xv n {x)v n ^ {x) 

h n -l,N J-oo A — X 2 

Let us prove that U n satisfies the linear equation 



* UnK : , n r ' M*)- (203) 



A(t/ n+ i — U n ) — r n+ i.N{U n +2 + U n+ i) — r n ^ N (U n + U n ^). (204) 
From (12031) we deduce that 

2 f°° 

\{U n+1 - U n ) = / xv n+1 (x)v n (x)dn(x) 



h 



n,N 



oo 

oo 



h n -i. 



N 



xv n (x)v n _i(x)&n(x) 



oo 



2 f°° x 3 v n+1 (x)v n (x) 

K,N J-oo X- X 2 

2 f°° x 3 v n (x)v n - 1 (x) 



h n -i,N 

Using x^v n (x) = LPv n {x) for j = 1, 2 we find that 



|°° iMM^). (205) 



xv n+1 (x)v n (x)dfj,(x) = h n+ltN , (206) 



and 

"™ x 3 w n+ i(x)t; n (x) 



d/x(x) = 



'CO 



d/i(x). 

A — x z 

(207) 

Substituting these identities in f !205j) and taking into account ( I203j) we con- 
clude that (12041 holds. 

It is now easy to prove that 

r„ lJV (£4 + £/n-i)(£4 + U n+1 ) = \{U 2 n - 1). (208) 



39 



Indeed, the linear identity (I204p implies 

r n+1>N (U n+2 +U n+1 )(U n+1 +U n )-r njN (U n+1 +U n )(U n +U n „i) = \(U* +1 -U*), 

(209) 

and therefore the expression 

\U 2 n - r niN {U n+1 + U n ){U n + U n -!) (210) 
is independent of n. Since r A r = and U = 1 the identity (12081) follows. 
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